Impact of the cracks lost in the imaging process on computing linear elastic properties from 3D microtomographic images of Berea sandstone

With the current developments in imaging/computational techniques and resources, computational rock physics has been emerging as a new field of study. Properties of rocks are examined by carrying out extensive numerical simulations on rocks that have been digitized using high-resolution X-ray CT scans. The ultimate goal of computational rock physics is to supplement the traditional laboratory measurements, which are time consuming, with faster numerical simulations that allow the parameter space to be explored more thoroughly. We applied the finite-element method to compute the static effective elastic properties from 3D microtomographic images of Berea sandstone saturated with different fluids. From the computations, we found discrepancies between the numerical results and the laboratory measurements. The reason for such a problem is the loss of small features, such as fine cracks and micropores, in the digitized matrix during the imaging and aggregation process. We used a hybrid approach, combining the numerical computation and the effective media theories — the differential effective medium model and the Kuster-Toksöz model— to deduce the lost cracks by a very fast simulated annealing method. We analyzed the sensitivity of the inverted results — the distributions of crack aspect ratios and concentrations — to the clay content. We found that the inverted crack distribution is not so sensitive to clay content. Compared with the effect of cracks on the computed effective elastic properties, clay has only a secondary effect. Our approach can recover the lost cracks and is capable of predicting the effective elastic properties of the rocks from the microtomographic images for different fluid saturations. Compared with the traditional inversion schemes, based only on the effective media theories, this hybrid scheme has the advantage of utilizing the complex microstructures that are resolved in the imaging process, and it helps define the inversion space for crack distribution.